Effect of the Coulomb interaction in A(d,p) fragmentation 
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In the framework of Glauber-Sitenko model we calculate contribution of Coulomb interaction 
in cross-section of A(d,p) reaction at high energy and zero angle. It is demonstrated that such 
effectsignificantly increases the differential cross section only at peak, where the proton momentum 
p is near half of the deuteron momentum pd in lab. frame, p ~ \pd- The Coulomb interaction 
do not change the results in the high momentum region, where quark effects should be taken into 
account. 

PACS numbers: 



I. INTRODUCTION 



For several decades considerable efforts have been done 
to investigate structure of the deuteron in a wide inter- 
nuclear region, from that where description is given in 
terms of nucleons and mesons to that where quarks and 
gluons should be explicitly used for the deuteron descrip- 
tion. Reactions with "hadron probe" , due to its high lu- 
minosity, play important role in such studies. Here we 
mention only the A(d, p) break-up at zero proton angle 
0, 1 1 1 E S 0, 1 M M M and the elastic pd- 
scattering at 180° in cm. frame 13 , [Til, Ell ■ For both re- 



actions detailed data on differential cross section and po- 
larization obscrvables, tensor analyzing power, T20, and 
polarization transfer, kq, were obtained. 

This data show significant deviation of the observ- 
ables (cross-section, T20 and kq) from theoretical calcu- 
lations, which cannot be removed by multiple scatter- 
ing effects without modification of the short-range part 
of the deuteron wave function [l6|. Such modification 
comes from Pauli principl on a level of constituent quarks 
[13, [3 . As a result, the short range part of the deuteron 
wave function includes, apart from the pn component, 
new components, NN*, N*N and N*N*, which cannot be 
reduced to the pn configuration. Because the main con- 
tribution of such states comes from the lowest resonances 
which have negative parity, the modification generates ef- 
fective P-wave in the deuteron which drastically changes 
behaviour of the observables at high momentum region. 

In the framework of such approach a good description 
of all break-up data was obtained [l6[. Nevertheless one 
important problem was not yet discussed. The data were 
measured at zero proton angle, 6 p = 0°, and Coulomb 
interaction might, in principle, give a sizable contribution 
in the observables. In the present analysis, we account 
for Coulomb interaction by adopting the Ahkiezer and 
Sitenko approach [lj[ [28| to inclusive 12 C(d,p) break-up 



and find that it significantly increases the cross section 
only at peak when the proton momentum p is near half of 
the deuteron momentum pd in the lab. frame, p ~ \pd- 
The paper is organized as follows. In Sect. [IT] we for- 
mulate the general formalism. We get expression for the 
Coulomb corrected amplitude (Sect. lH~A"|) . give remarks 
to strong interacting part of the amplitude (Sect. Ill Bp 
and obtain the general expression for the invariant dif- 
ferential cross section of the (d, p) break-up at zero proton 
angle (Sect. Ill C|) . The procedure of numerical calcula- 
tions is discussed in Sect. IIIII In Sect. IIVI we compare 
our calculations with experiment and give a summary. 



II. FORMALISM 

A. Coulomb corrected amplitude 

We begin by considering general formalism for the 
Coulomb correction to cross-section of the deuteron 
break-up A(d,p) at the proton longitudinal momentum 
P3 ~ \v<i and small transverse momentum p±. In this 
case it is enough to take into account the elastic contri- 
bution only. 

Neglecting the Coulomb interaction the amplitude for 
the elastic break-up reads 



^ Str (P-L,P3,Cu) = f^ Jd*B 

x / d 3 r?Pl(r)Mr) 



iXstr(b p ,b„; 



(1) 
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= f* / d 2 Be^ B J d 3 W4(r)Mr)x 

x [r„(b„) + T p (b p ) - r„(b n )r p (b p )] = 

= F n ^-( P± , P3 , Q ± ) + F p ^-(p ± ,p 3 , Q x ). 
-F„f-(p ± , P3 ,Q±), 



where ipo ( r ) and vb^ (r) are the wave functions of the 
deuteron and the final proton-neutron system and T p (h p ) 
and r„(b„) arc the profile functions for the proton and 
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the neutron; B = j(b p + b„) and = b p — b n , with 
bp and b„ are impact parameters for the proton and 
the neutron; = (0,0, pa) is the deuteron momentum, 
Q± is the momentum transferred to the final proton- 
neutron system. The definition of the relative momentum 
between the final proton and neutron, k, needs a special 
comment. In the non-relativistic case it is defined as k = 
(p_L — |Q_i_,P3)- For the relativistic deuteron, pd 3> ma, 
we define it by boosting along z axis to a frame where the 
total longitudinal momentum of the two-nucleon system 
is zero, 



p* 3 + n|=0, p!=P_l, n^=n ± , (2) 

k± = |(p_L - nj.) = p± |Qx, h = |(P3 - "aM 3 ) 

where p* and n* are momenta of the final proton and 
neutron in the new frame; p and n arc the same momenta 
in the lab. frame. Of course, such definition can be 
accepted only near the region p 3 ~ \pd- Assuming that 
the transverse motion is non-relativistic one simply gets 



v^ 7 



E d 



P, 



where 



I>5 



-p d +p. (4) 



The profile function is given by 
1 



r(b) 



2nip 



N 



d 2 le- Ub f N (q) 



(1 - iPN)crN itf/pl 

Air/3% 



where 



f N (l) 



(i + pn)pn&n c -±i3 2 M i 2 

47T 



(5) 



(6) 



with fUjf is the slope parameter for the nucleon-nuclcus 
scattering, a N is the total cross section and pm is the 
ratio of real to imaginary parts of the amplitude. Laler 
on we will use a natural approximation a p = <r n = a, 

f3 p = f3„ , = P and pn = 0. 



Following the prescription of Ref. [21( the Coulomb in- 
teraction is included by adding the Coulomb, Xc(b p ), and 
screening Coulomb, Xscr, phase shifts to strong interact- 
ing phase function x s tr(b p , b„), i.e. 

r(b p ,b„) = 1 - e %tr(b p> b n ) 



I _ gi[Xstr(bp,b n )+X B cr+Xc(6p)] 



(7) 



= r(b p ,b„) + e^< 

x [i-r(b p ,b„)], 

where the Coulomb and screening Coulomb phase shifts 
are given by 

., . 2Za 4itZu 
Xc{bp) = mp 3 H lno. 



T c (b')]nb'b'db' 



ib p 
2Za 



i> 



ln2p 3 i? scr , 



T c (b')b'db'+ 

Jo 

= Xo + Xc(b P ), (8) 
(9) 



respectively Here T c (b) = J p c (r)dz is the thickness 
function corresponding to the nucleon charge distribu- 
tion, p c {r), normalized by J p c (r)d 3 r — 1; v p is the pro- 
ton velocity, R scr is atomic screening radius, Z is the 
atomic number and a w 1/137 is the fine structure con- 
stant. 

After that the Coulomb amplitude reads 



F C ( PJ ., P3 ,Q±) 



Wd 
2tt 



d 2 Be *Q ± B / d 3 rx 



(10) 



x f k (r)^ (r) | 

x (i-r p -r„ + r p r„) = 

dis p n 1 pn 

Now let us consider contribution of different terms of 
this expression. 

• Coulomb dissociation. One simply gets 



^ c is =f^(iQ,k)x 



d 2 b p e tbpQ 



where 



G(iQ,k) = J d 3 re^^i(r)Mi 



(11) 



(12) 



is the transition form factor. To perform integra- 
tion over d 2 br, it is useful to introduce a point charge 
phase shift [211 ] 



X P t( 
where 



2Za 



lnp 3 b p = xo + X P t(bp), (13) 



X P t(M 



2Za 



In b D 



and rewrite identically [e 



p'Xpt 



] + [e 



arrives at 

FS* =G (|Q, k) [^(Q) + AF(Q)} ee 
=G(±Q,k).F c (Q), 
where 

2Zap d 



Finally one 
(14) 



^pt(Q) 



vQ 2 



(15) 



is the scattering amplitude on Coulomb potential of 



a point charge (<p c = — — 
C = 0.577215 is Euler constant) and 



, where 



poo 

= ipd / dbpbpJo(Qb p ) 



i' 



where Jq(x) is the Besscl function. 



JXc{b P 



(16) 
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• Proton scattering affected by Coulomb interaction. 



^ C = ^G(-iQ,k)C c (Q), 



(17) 



where C C (Q) = dqqe-^+Q^ I Q {^Qq)T c {q)- 

• Neutron scattering affected by Coulomb interaction. 

F: = J ^G(iQ-l,k)/ n (0^ c (|Qx-lx|). (18) 

• n — p rescattering affected by Coulomb interaction. 



[ d 2 ld 2 l'G(-±Q + l',k) x 
rj^pn J 



(2tt) 2 

X^ c (|Q-l-l'|)/ p (0/n(O- 



(19) 



For gaussian amplitudes f p (l) and f n (l') one can 
integrate over the angles 



dq / dq'qq'G{±tf,k) X 
JO 



xT c (q)I (-y 2 qq')l (^ 2 qQ) X 

x 



(20) 



where Io(x) = e~^Io(x). Deriving (12"0)) we used 
the fact that the form factor G(Q', k) does not de- 
pend on the Q-direction when the vector k is di- 
rected over z axis. We have also assumed that the 
nucleon-nuclei amplitudes are the same for the pro- 
ton and the neutron. 



Finally one gets 



™tr. , i Xs0 r (pc _pc_ F c,pc\ 
± ~ c V dis ± p ± n ' ± npj 



(21) 



As the result of infinite range of the Coulomb inter- 
action, the amplitude (fl~5|) is diverged at Q — > 0. Due 
to the deuteron bound energy the neutron loses a part 
of its longitudinal momentum Qn at the limit when the 
transverse momentum Q± — > [19|. We estimate it to 
be 



Qo = 



m 2 n + (Pd-p) 2 -{E d -E p f 
2(pd-p) 



(22) 



To take into account this effect, one has to change Q 2 — > 
Q 2 + Qa m the denominator of the amplitude ((T5|) . The 
momentum Q is minimal at p = ^pd (see Figure 1) 
what means that at this kinematical region the Coulomb 
interaction should be maximal. 



constituents (the proton and the neutron) suffer elastic 
scattering only. Its square is proportional to the so-called 
"disintegration cross section" . To calculate the inclu- 
sive cross section, one has to add two contributions, the 
disintegration and absorption cross sections [13, HH . In 
the later the neutron suffers inelastic collisions, but the 
proton keeps elastic scattering only. It is the core of 
Bertocchi-Treleani model. 

In Ref. [l6| Bertocchi-Treleani model (23| was modified 
in the following way: 

• the deuteron wave function was considered in the 
framework of "minimal relativization prescription" 
with dynamics in the infinite momentum frame [24|, 
H; 

• it takes into account Pauli principle at the con- 
stituent quark level. 

According to Resonating Group Method (RGM) Pauli 
principle at the level of constituent quarks modifies the 
deuteron wave function at short distances by 0, EH 

^(1, 2, 6) = 1<M1, 2, 3)<M4, 5, 6) X (r), (23) 

where A is the antisymmetrizer for quarks from different 
three quark (3g) bags, 4>n are nucleon 3q clusters, x(r) is 
the RGM distribution function and r stands for the rela- 
tive coordinate between two 3g bags. Due to the presence 
of the antisymmetrizer in (|23|) , the deuteron wave func- 
tion, being decomposed into 3q x 3q clusters, includes, 
apart from the standard pn component, nontrivial NN* 
components. Most of the isobars N* have negative parity 
and thus they generate effective P-wave in the deuteron. 

Following Refs. [Tt], one can choose \( r ) as a con ~ 
ventional two nucleon deuteron wave function renormal- 
ized by the condition, formulated in Ref. [26j . 
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B. Strong interaction part of the amplitude 

As was mentioned before, the strong interacting ampli- 
tude (Q]) corresponds to the situation when the deuteron 



FIG. 1: Dependence of the Qo momentum (|22p on the proton 
longitudinal momentum p for the deuteron momentum pd = 
9.1 GeV/c. The solid line is for the experimental masses of 
the proton and the neutron and the dashed line is for m p = 
m n = hmd- 
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The cross section 



The differential cross section is given by 

P k =^sJ^\nP„PM\ 2 , (24) 

where p and n are the proton and neutron momenta. In 
the case, when p± = the transverse momentum Q = 
and the integral over the angle becomes trivial. Finally 
one arrives at the following expression for the invariant 
differential cross section 



In numerical calculations we use parametcrizations for 
the phase shift Xc(b) © and the real and imaginary parts 
of the potential correction AJ r c (Q) defined in Eq. (fTo)) 
which are given in Appendix. 

All calculations arc done with the dcutcron wave func- 



tion for Nijm-I potential [27j. The 5 and D wave com- 
ponents of the wave function were approximated by sum 
of gaussians 



N 



u(r) = r Aj 



w(r) 



-0ir 2 



(28) 



i=l 



<Pp (2tt)2 J 



dnn\F(0 ± ,p 3 ,n)\ , (25) 



where E* is the proton energy in the deuteron rest frame. 

Due to oscillation factor e IXsor in (|2"Tj) . the strong in- 
teracting and Coulomb parts of the amplitude do not 
interfere in (l25l). 



III. NUMERICAL CALCULATIONS 



For the gaussian wave function (|28|) one cannot construct 
wave function for the unbound pn system which fulfill the 
conditions of orthogonality and complitnes. Similar to 
[20| one can construct only a function which is orthogonal 
to the deuteron wave function 



V> k (r) = e lkr - (2^) 3 / 2 Mr)Mk)/N s , (29) 

where Vs( r ) = ^y/2 r u ( r )i N s is probability for the S- 
wave in the deuteron and <p s (k) is the Fourier transform 
of the deuteron S'-wave function, 



We use the nuclear charge density which corresponds 
to the harmonic oscillator well 



Pc(r) 



IT 2 Z<Xq 



1 + (Z - 2) 3^ 



(26) 



>(*) 



FT N 

— Y 



1 (2a,) 3/2 



-fc 2 /(4a.) 



Finally one gets the following expression for the form 
factor 



where (r 2 



Z (2 4, 



For 12 C the parameter gq = 



1.60 fm 12 111 . The thickness function 



T c (b) 



ttZciq 



Z -2 



OUn 



i.2 / 2 

-6 /o 



(27) 



The Coulomb phase shift Xc{b p ) is displayed at Fig- 
ure 2. 



0.2 



G(iQ,k) = (2tt) 

N 



3/2 



(l|Q-k|)- 



T 3/2 w) 



AiAj 



i-.r- 



■ exp 



Q 2 



16(a, 



(30) 



In numerical calculations we take for the total p 12 G cross 
section its experimental value a =340 mb. The slop 
parameter was calculated in the framework of Glauber- 
Sitenko model to be {? = 69.3 (GeV/c) 2 . 
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FIG. 2: The Coulomb phase shift x c (b p ) (the solid line) and 
the phase shift for a point charge Xpt(b P ) (the dashed line) 
for 12 C. The dashed-point line is for the approximation of the 
phase shift x a p.(M- 



IV. COMPARISON WITH EXPERIMENT AND 
SUMMARY 



In Figure 3 the results of calculations are compared 
with experimental cross section data. We also com- 
pare our results with quasi-impulsc approximation, i.e. 

ss <T n I cl' ! /'d( < 7)| 2 ; where q is the proton momentum in 
the deuteron rest frame and <j™ c is the inelastic neutron- 
carbon cross section[29]. One sees that the Coulomb in- 
teraction strongly inrcases the cross section at the peak 
near q ~ (which corresponds to p ~ \pd)- Out of the 
peak such effect sharply decreases and becomes negligibly 
small at q >100 McV/c. So it cannot affect the region, 
q >200 MeV/c, where quark effects are assumed to be 
significant. 
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FIG. 3: The invariant cross section of the 0° inclusive 
12 C(d,p) break-up at pd = 9.1 GeV/c plotted versus the pro- 
ton momentum in the deuteron rest frame. The experimen- 
tal points are from [J 0| • The dash-dotted curve is for the 
quasi-impulse approximation (see text), the dashed curve is 
for multiple scattering + Pauli principle at quark level (PQL) , 
the full curve is for multiple scattering + PQL corrected by 
Coulomb interaction. 
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APPENDIX A 



In numerical calculations we use the following param- 
eterizations for the phase shift Xc(b) 



X ap .(6) = — In | 
v p 



A, 



1 + A 2 b 2 



(Al) 



(where A 1 =2.2416 fm 2 , A 2 =0.34 fm" 2 and b is defined 
in fm) and the real and imaginary parts of the potential 
correction AjF c (Q) 



KeA^ c (Q) = - 
9fmA^ c (Q) = - 



where ip r = 
0.035 GeV 2 , B r 



2Zap d 
v 

2Zap d 

v 



2Zap d 



COS if r 

A r + Q 2 
Aj 



+ C r Qe 



Qe 



(A2) 



C + lnf=r±2 



2p 3 



and A r = 



0.086 GeV, C r = -12.410 GeV" 



0.125 GeV^, A; 



0.264, Bi 



0.799 GeV" 



Qf = 0.0242 GeV , C, = -1.112 GeV" 
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